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ABSTRACT 

We present a new near-infrared photometric system for detection of water 
ice and methane ice in the solar system. The system consists of two medium- 
band filters in the K-hand region of the near-infrared, which are sensitive to 
water ice and methane ice, plus continuum observations in the J-band and Y- 
band. The primary purpose of this system is to distinguish between three basic 
types of Kuiper Belt Objects (KBOs) — those rich in water ice, those rich in 
methane ice, and those with little absorbance. In this work, we present proof- 
of-concept observations of 51 KBOs using our filter system, 21 of which have 
never been observed in the near-IR spectroscopically. We show that our custom 
photometric system is consistent with previous spectroscopic observations while 
reducing telescope observing time by a factor of ~ 3. We use our filters to identify 
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Haumea collisional family members, which are thought to be colhsional remnants 
of a much larger body and are characterized by large fractions of water ice on 
their surfaces. We add 2009 YE7 to the Haumea collisional family based on our 
water ice band observations (J — H2O = —1.03 ± 0.27) which indicate a high 
amount of water ice absorption, our calculated proper orbital elements, and the 
neutral optical colors we measured, V — R = 0.38 ±0.04 , which are all consistent 
with the rest of the Haumea family. We identify several objects dynamically 
similar to Haumea as being distinct from the Haumea family as they do not have 
water ice on their surfaces. In addition, we find that only the largest KBOs 
have methane ice, and we find that Haumea itself has significantly less water ice 
absorption than the smaller Haumea family members. We find no evidence for 
other families in the Kuiper Belt. 

Subject headings: comets: general — Kuiper Belt — solar system: formation 



1. Introduction 

Detecting ices on small outer solar system bodies is of utmost importance because the 
presence of the most volatile ices, such as methane ice, may indicate primordial surfaces. 
The most massive of the known Kuiper Belt Objects (KBOs) such as Eris and Pluto, cannot 
retain surface ic es for the age of the solar system even under moderate heating, temperatures 



of > 70 Kelvin f lSchaller fc Brownll2007bl ). Thus, an inventory of volatile ices on the KBOs 
has the potential to provide a useful metric of historical heating. In the absence of planet- 
wide resurfacing, those KBOs that originally had ices but passed too close to the sun would 
have lost their ices and those which stayed distant could retain their ices. Tying the presence 
of surface ices to KBO orbital classes is one of the few direct observables which can shed 
light on the principle dynamical mechanisms involved in shaping the early Kuiper Belt. 



^ Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United States), the Science and Technology 
Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and Ministerio de 
Ciencia, Tecnologia e Innovacion Productiva (Argentina). 

^This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas 
Observatory, Chile. 
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Most ices of interest are bland or nearly so at visible wavelengths, but many show 
prominent absorption features at near-infrared wavelengths. Although many ices show strong 
features in the mid-infrared, at these wavelengths remote sensing via ground-based or space 
based telescopes is severely hampered by thermal background radiation, making detection 
of ices difficult. Thus, nearly all surveys of KBOs to date that have positively identified 
ices have used near-infrared data to do so. The only major sur face compon ent which is not 
detectable in the near-infrared is the so-called ultra-red matter (IjewittI [20021 ). which appears 
in the visible and has not been definitively shown to be correlated with near-infrared colors. 

Many researchers have collected near-infrared spectra of the brightest KBOs. Of the 
~ 30 spectra of these bodies, almost all can be categorized into three types: water ice dom- 
inated, methane ice dominated and neutral. For example, the five brightest KBOs (plus 
Pluto and C haron) can be div i ded into either wa t er ice dominated (Charon, Haumea, Orcus 
and Q uaoar: | Buie et all JiQStI): iBrown fc Calvin! J2000h: Ijewitt fc Luul fl2004j ): iTruiillo et al. 



(120051 . l2007al ): 



Pinilla-Alonso et al. 



(120091): 



Eraser fc Brown 



(2009)) or methane ice domi- 



nated ( Pluto, Makernake, E ris; Cruikshank et al. ( 19761 ): Brown et al. ( 2007a . 2005 ): Dumas et al 



( 120071 ): iMerlin et al.l ( l2009l )). Many of the fainter bodies have been shown to be featureless. 
There are some bodies that fall outside this broad taxonomy and have combinations of 
materials such as T riton, which show significant absorption of combinations of substances 
( lOuirico et al.l 119991 ). Such bodies may not be adequately described by a simple set of pho- 
tometric information. We note however, that our survey goal is not to replace the detailed 
information that spectroscopy can yield, but to be able to identify population-wide trends 
and KBOs with unusual surfaces for follow-up spectroscopy. Although most bodies can be 
grossly characterized as either water or methane types, it is still important to note that many 
ices remain extremely difficult to detect in the near- infrared such as carbon monoxide (CO), 
carbon dioxide (CO2), and nitrogen (N2) due the characteristically narrow and weak absorp- 
tion lines of these ices and the extreme faintness of the KBOs. Even for Pluto, which is over a 
factor 10 brighter than the ne xt brightest KBO, detection of CO, CO2 and N2 is challenging 
with modern instrumentation ( lOlkin et al.ll2007l ). Eor this reason, tracing methane ice in the 
Kuiper Belt is extremely important as it is one of the few volatiles that is easily detectable 
and very sensitive to thermal history. Water ice acts as an important comparison ice, as it 
is not volatile beyond ~ 5 AU in the current solar environment, although like methane it 
can be destroyed by galactic cosmic rays and solar UV photons, likely two of the dominant 
long-term weathering processes in the Kuiper Belt (jCooper et al.l 120031 ). 



There have been several photometric programs, with varying levels of s uccess. Numerous 



attempts at JHK photometry have been published (see review work by lEulchignoni et al. 



20081 . and references therein), but unfortunately, these bandpasses are only serendipitously 
correlated with methane and water ice transitions and cannot readily discriminate between 



-4- 



the tw o. Probably the most successful set of observations to date are those of lSnodgrass et al. 



( 120101 ). who surveyed Haumea family members in search of water ice in the if-band. Their 
use of the medium bandpass Hg filter (1.52 - 1.63 /im) which probes the solid-state water 
ice region in KBOs allowed the identification of many Haumea family members using their 
photometry and references to literature. Although the Hs filter is sensitive to water ice, it 
is much less sensitive to methane ice, which has peak absorption between 1.60 - 1.85 /im in 
the if-band. 

In this work, we introduce a simple custom-bandpass three filter system which allows us 
to discern between KBOs with water ice, methane ice and relatively bland surfaces. We have 
undertaken proof-of-concept observations which demonstrate that these filters are indeed 
consistent with spectroscopic observations yet typically require 1/3 of the telescope time of 
spectroscopic works. In a ddition, we have observed all of the Haumea family candidates 



brighter than = 22 from iRagozzine fc Brownl (120071 ) and we identify which members have 
surfaces rich in water ice. 



2. Experimental Design 

2.1. Custom Filter Bandpasses 

Our filter set was chosen by simulating synthetic spectra of water ice, methane ice 
and neutral KBO surfaces as observed through the Earth's atmosphere using data from 
the Gemini/NIRI Integration Time Calculator to account for instrumental transmission and 
thermal emission. For the methane ice filter, regions in the K-hsnd were favored over the 
if-band for 4 reasons: (1) the iJ-band contains many more strong telluric OH emission lines 
than the i^-band, which are variable on ~ 5 minute timescales and may affect photometry, 
(2) transitions for methane and water in the i^-band are deeper than in the i7-band and are 
more widely separated from each other, (3) a wider bandpass without strong atmospheric 
extinction is available in the i^'-band and (4) the Gemini North telescope, used for the 
faintest bodies in this program, is optimized for near-infrared observations and thus has 
low emissivity and low thermal contamination in the i^'-band. The possibility of double- 
bandpass "notch" filters including ice transitions in both the if-band and the /T-band while 
excluding the sky between were also considered, but estimated signal gain was not sufficient 
to warrant the inclusion of the i7-band telluric OH emission and the loss of throughput for 
the double-bandpass design. 

In order to estimate the amount of absorption seen in the filter, continuum observations 
(i.e. observations out of absorption bands) are also needed as a reference point. Continuum 
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observations must be in regions that have minimal (1) telluric absorption, (2) atmospheric 
emission and (3) ice absorption from the target body. As seen in Figure [H the optimum 
region that satisfies all 3 criteria is in the J-band, where water, methane, and atmospheric 
contamination are minimal. There are no regions in either the if-band or fC-band which are 
devoid of atmospheric, water and methane ice absorption. We also collected F-band data on 



many bodies as a mitigation against the presence of so called "ultra-red" material (jjewitt 



2OO2I ). which can extend throughout the visible and contaminate the shortest wavelengths of 



the J-band. After final analysis, robust results were obtained without the use of the F-band 
filter data, but in a few cases, such as Quaoar, which is extremely red in the visible, it was 
needed (see Section [5.2.ip . Final filter bandpasses are listed in Table [TJ 



2.2. Filter Manufacturing and Installation 

We installed our filters into two telescopes — The 8.1 m mirror diameter Frederick C. 
Gillett Gemini telescope located atop Mauna Kea, Hawaii and the 6.5 m mirror diameter 
Magellan Baade Telescope, at Las Campanas, Chile. Our filters were manufactured by Barr 
Associates, Inc. in Massachusetts, USA, who implemented the bandpasses to our optical 
specifications. Filters with identical bandpasses were manufactured in the same coating 
batch to ensure consistent quality — only the geometric shape of the filter substrate differed 
between Gemini (60 mm diameter, round, 5 mm thick) and Magellan (32mm across, square, 
4 mm thick) . The filter efficiencies are presented in Figure [2l Filters were i nstalled in late 



2008 in the Gemini North Near InfraRed Imager and Spectrometer (NIRI, iHodapp et al. 



2003f) and in the M agellan Baade Persson's Auxiliary Nasmyth Infrared Camera (PANIC, 



Martini et al.ll2004j ) dewars. Both Gemini/NIRI and Magellan/PANIC are cryogenic instru- 
ments, so filter installations are performed very rarely, due to the minimum of 2 weeks warm 
up and cool down cycle time required to service the instruments. An extra set of Gem- 
ini/NIRI filters was manufactured to provide contingency in case unexpected manufacturing 
defects or installation problems were encountered. 



2.3. Filter Commissioning and Transformation to Common Bandpasses 

On-sky tests of the filters were performed shortly after installation to ensure acceptable 
quality of the filters. These tests involved imaging a UKIRT JHK standard star to estimate 
throughput and zeropoint and to ensure that image quality was consistent with other stan- 
dard filters. A full characterization with additional targets has since been performed with 
a large number of UKIRT standards observed over several semesters at Gemini North. We 
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present our custom color metrics versus standard color metrics in Equations [T] and [2l which 
are derived from the data displayed in Figure |3l 



It is important to note that while these transformations are applicable to UKIRT stan- 
dard stars which generally have few transitions in the K-hand, they are not applicable to 
KBOs or other bodies with icy surfaces, as the entire purpose of our experiment is to probe 
the absorption regions of KBO ices. 

Since the Gemini North NIRI instrument is also a spectrometer, spectra were collected 
of the filters at cryogenic temperatures indicating that the overall throughput was consistent 
with that measured in the Barr laboratory (Figure [2]) — thus no degradation of filter per- 
formance occurred during installation or cooling to cryogenic temperatures. This was not 
possible to do with Magellan PANIC since it is an imager only, however, the Magellan filters 
were manufactured using the same techniques and in the same batch as the NIRI filters, so 
we expect that the filters are identical in their physical properties. Results of a few science 
targets observed with both instruments showed consistency. 



Near-infrared observations of the brightest targets were performed at Magellan on four 
nights spanning two observing runs on UT dates 18-19 October 2008 and 03-04 April 2009. 
Weather conditions were photometric for both runs with typical delivered image quality 
through our custom filters of 0.3 - 0.5 arcsec. Gemini observations were conducted in the 
Gemini North Queue under program IDs GN-2008A-Q-4, GN-2008B-Q-40, GN-2009A-Q-6, 
and GN-2009B-Q-30. Observing conditions were set at 70 percentile or better for image 
quality (i.e. i^-band image quality of 0.55 arcsec or better at zenith) and 50 percentile for 
cloud cover (i.e. photometric). This typically led to image quality similar to that seen at 
Magellan, from 0.3 - 0.5 arcsec delivered on the detector. A more complete summary of the 
objects and observational circumstances of the observations appears in Table |2l 

At Magellan, the telescope was tracked at the sidereal rate, while at Gemini, the tele- 
scope was tracked at the non-sidereal rate of the object. These two methods were used for 
ease of execution at the telescope — at Magellan any non-sidereal tracking rates must be 
updated upon slew by hand while at Gemini the telescope can automatically compute mov- 
ing object position and tracking rates from orbital elements. Since we were limited to ~ 60 



J-H2O = 0.991(J 
J-CH4 = 1.104(J 



K) -0.313 
K) - 0.367 



(1) 
(2) 



3. Near-Infrared Observational Program 
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second exposures or shorter due to the thermal sky background, our objects, which typically 
have ~ 3 arcsec / hour apparent motion, moved less than 1/5 of a seeing disk during each 
image, resulting in image quality being identical for both the object and background stars 
no matter which guiding technique was employed. 

While observing a science target, filters were cycled with repeated visits to the J- band 
to allow for identification and correction of any lightcurve effects. A typical observational 
sequence for a faint target at Gemini was 3 images taken with 30 second exposure time (3 
X 30s) J, 6 X 60s H2O, 6 x 60s CH4 repeated for 1.25 hours while offsetting the telescope 
between each exposure. In addition, a short sequence of 3 x 30s Y was included in the 
beginning of the sequence as an additional constraint on any large scale color trends that 
might be visible in the 1/xm wavelength region. 

Primary color calibration was done by observing a G2V solar analog star immediately 
before or after every science target at an airmass similar to the science observations. The aim 
of collecting these data was to estimate any telluric extinction or emission in the custom filter 
regions, which was particularly important in the H2O filter which contains a known region 
of atmospheric attenuation, which could result in a ~ 3% photometric variation (~ 0.03 
magnitudes) if left uncorrected. These techniques are similar to what is routinely done for 
spectroscopic telluric correction to provide valid spectral information in the 2//m wavelength 
region. This provides for the simplest possible measure of object reflectance — science target 
colors arc compared to the G2V analog taken at similar airmass to immediately and directly 
calculate reflectance. 

In addition to the G2V solar analog, a UKIRT standard star was observed every night 
science observations were conducted. This provided information about any possible long 
term degradation of the fllters which could affect our results if it occurred. No photometric 
change in filter sensitivity of color was observed over the course of our observations. 

4. Data Reduction 

Basic data reduction was performed in a standard manner for near-infrared imaging and 
in a similar manner at both telescopes. First, daytime flat flelds were examined to identify 
bad pixels (i.e. those pixels that did not respond linearly with a change of flat fleld level). 
These pixels were replaced with a median of nearby pixels on all science and calibration 
images. Flat fields at Gemini were constructed from daytime observations of a calibration 
source, while at Magellan sky frames were used. The science fields were divided by these 
fiats to correct for the pixel-to-pixel variations seen in the detector. At Magellan, the pixel- 
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to-pixel variations were very minor as the detector was generally cosmetically pristine for 
a near-infrared array. Two deep images were then created for each target, consisting of all 
science frames taken of a particular object (1) digitally shifted to track the object at the 
non-sidereal rate of the object of interest, and (2) digitally shifted to track at the sidereal rate 
to create a deep image of the star field. This allowed for the rejection of science observations 
where the science object encountered background stars during the course of the observations. 
Finally, synthetic aperture photometry was performed on all of the science and calibration 
targets. 

Photometric calibration was performed on a nightly basis by identifying all calibration 
stars taken on a given night and plotting the stellar brightness as a function of airmass. Any 
deviant calibration targe ts were rejected as b eing stars misidentified as G2V in the Simbad 
Astronomical Database (jWenger et al.ll2000l ). with spectral typing derived from the Mich i- 
gan Catalogues of Two-Dimensional Spectral Types for HD stars (IHouk fc Cowleyl Il975l ). 
Remaining stars were used to create a telluric correction curve as a function of airmass. 
In general, these observations remained consistent on a nightly basis. These photometric 
calibrations were applied to the science observations to produce a relative reflectance mea- 
surement for each object. 



5. Results 

The results of our work are presented in Figure |H The basic goal of this work is to 
demonstrate that our custom filter bandpass can distinguish between the three basic types 
of bodies observed in the Kuiper Belt: water ice, methane ice and neutral reflectance. This 
trend is readily apparent in Figure H] as the vast majority of objects fall into the neutral 
category in the upper right. Most bodies with ice detection have water ice and occupy the 
left part of Figure HI The two known bodies with methane ice, Makemake and Eris are 
shown in the bottom right of the figure, clearly separated from all other objects. 

In this proof-of-concept sample, we have studied the brightest KBOs for two reasons: 
(1) these bodies require less telescope time to study and (2) many of these bodies have 
spectroscopy already collected by other researchers which provides corroborative evidence 
that our survey is performing as expected. We believe that the current sample adequately 
demonstrates that our custom photometric filter system does indeed distinguish between 
these basic types of bodies as the three basic surface types are clearly discernible by eye. 
In future works, we will extend our survey to fainter bodies requiring more telescope time 
— bodies that to date have been inaccessible to spectroscopic studies. Even in our current 
sample, about 40% of our targets have no known spectra collected, primarily due to their 
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faintness. 

Table E] provides verification data for tliis work. It lists the color quantities computed 
for each object and the surface type based on this color quantity. In addition, it also lists 
references to published works that provide spectroscopic verification of our basic surface 
identification. Discrepancies between the two methods occur only for bodies with low ab- 
sorption, as discussed in§ 15. 2[ All of the bodies with strong water ice features show excellent 
correlation with published spectra. The identification of water ice among the Haumea fam- 
ily members, including recently discovered object 2009 YE7, as well as Ixion and Orcus, are 
discussed in section § 15. 1[ For the Haumea family, we find that the faintest bodies system- 
atically have the strongest transitions, as discussed in § 16. 1[ In all cases, our photometric 
study was done in far less time than spectroscopic identifications, as discussed in further 
detail in § O 



5.1. Family Members 

5.1.1. Haumea Family 



The Haumea family members are readily apparent by examining Figure H] and cross- 
correlating icy bodies with orbital parameters. We present this cross-correlation in Table IH 
where we list all bodies with J — H2O < —0.3 (water ice bodies) or J — CH^ < —1.3 
(methane ice bodies). All of the Haumea family members are readily identifiable as having 
both water ice on their surfaces and having proper orbita l elements, i ~ 28°, a ~ 43.5 
AU and e ~ 0.1, as presented in iRagozzine fc Brownl (120071 ). We have observed all bodies 
brighter than V = 2 2 listed as possible Haumea family members on dynamical grounds from 



Ragozzine fc Brownl (120071 ) and report physical associations based on their high percentage 
of water ice, as summ arized in Table El These identifications match previous spectroscopi c 
works on the subject (IBrown et al.l l2007bl : ISchaller fc Brownl l2008l : ISnodgrass et al.ll2010l ). 
We also report the identification of a possible new family member, 2009 YE7, based on its 
unusually high water ice fraction and proper orbital elements which are similar to other 
Haumea members. We computed the proper orbital elements for 2009 YE7 using an average 
of 20 million yea rs of orbital motion computed by the MERCURY symplectic integrator 
Jchamberslligggi ). 
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5.1.2. 2009 YEt 



After determining that 2009 YE7 was a Haumea family member, we proceeded to collect 
visible photometry of it to determine if it was consistent with the other Haumea family 
members which are neutral to slightly blue in reflectance, indicating not only that large 
amounts of water ice are on their surfaces, but that the water ice is in a largely pristine 
form. 

Optical observations of 2009 YE7 were obtained at The Magellan Clay 6.5 meter tele- 
scope on UT March 21, 2010. The LDSS3 camera has a STA0500A 4064 x 4064 Charge- 
Coupled Device with 15/im pixels. The field of view is about 8.3 arcminutes in diameter with 
a scale of 0.189 arcseconds per pixel. Each image was reduced using dithered twilight flats 
and biases in a standard manner. Images were acquired through Sloan r' and i' filters 
while the telescope was auto-guiding at sidereal rates using nearby bright stars. Exposure 
times were between 250 and 300 seconds. Souther n Sloan standard stars were used to pho- 
tometrically calibrate the data (jSmith et al.ll2005l ). In order to more directly compare our 
results with previous works, the Sloan colors wer e converted to the J ohnso n- Morgan- Cousin s 
BVRI color system using transfer equations from lSmith et al.l (120021 ). as in ISheppardI ( I2OIOI ). 
These colors are presented in Table [61 

The optical colors of 2009 YE7 have a spectral slope of S { i' > q ') = 5.3 ± 3, as defined 



m iDoressoundiram et al.l ( 120081 ) and implemented in ISheppardI ( I2OIOI ). Such a spectral slope 
is consistent to the other Haumea family members, which have a weighted mean slope of 
S' = 1.9 ± 0.3, as computed from the values presented in ISnodgrass et al.l (l2010l ). As the 
visible co lors of 2009 YE ? are inconsistent with most other KBOs, which have a mean 
~ 20 (ISheppardI I2OIOI ). we consider the visible colors to be confirmation of the near- 



infrared photometric and dynamical identification of 2009 YE7 as a Haumea family member. 



5.1.3. Other Families 



We find two objects in our survey have much more water ice on their surface than the 
other non-Haumea family members, namely Orcus and Ixion. These two bodies have low 
levels of water ice compared to the Haumea family members, as found by our s urvey and cor 
roborated by spectroscopic studies which show possible detection of water ice ( iTrujillo et al. 
2OO5I : iBarkume et al.ll2008l ). These two bodies both occupy the 3:2 mean motion resonance 
with Neptune and have inclinations within 1 degree of one another. Although a shared origin 
might be invoked for these bodies, adding in results of previous spectroscopic works suggests 
otherwise. Spectroscopic detections of water ice show no apparent correlation with inclina- 
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tion among the 3:2 resonance, as (47171), (55638), (84922) and (208996) have inchnations of 
8.4, 16.3, 14.8, and 13.6 degrees, respectively (Tableland references therein). Thus, at the 
present time, we find no evidence for a family in the 3:2 Resonance, although clearly there 
are bodies with moderate water ice present. 

The 3:2 Neptune resonance is thought to have "swept" up many bodies shortly after 
the solar system formed due to Neptune's outward migration ( lMalhotralll993l ). Thus, even 
if there were a collisional family in the 3:2 resonance, there is likely to be many non-family 
members within the 3:2 resonance that have no water ice on their surfaces. Identifying a 
family among interlopers would be very difficult in the 3:2 resonance unless the collision 
happened very recently, of order ~ 100 kyr, over which time orbital elements can be sta- 
ble for some objects. Much longer resonance habitation or large collisional velocities can 
cause prof ound departure (and ejectiori) of bodies, due to the very fine structure in the 3:2 
resonance (ITiscareno fc Malhotral l2009l : iMorbidelli et al.l 120081 ) . 



It has been suggested that KBOs may be devoid or somewhat depleted of ices due 



to bo mbardment from Solar U V photons and cosmic rays ( iMoroz et al.l |2004J : ICooper et al. 



20031 ). These ionizing particles are expected to destroy the molecular bonds of simple ices, 
resulting in their depletion, and eventually yielding a dark, carbon covered surface. Indeed, 
many faint KBOs have been studied spectroscopically, and in this work photometrically, the 
majority of which are neutral in refiectance. However, the presence of both methane ice 
and water ice, although the minority for all KBOs, is pervasive throughout the dynamical 
classes in the Kuiper Belt. In our survey, we find a Scattered KBO with ice (Eris), Resonant 
bodies (Ixion and Orcus) and Classical (Makemake) bodies with ices as shown in Table HI 
Among previously published spectroscopic works the same pervasive nature of water ice is 
apparent from Table |3]and references therein, with Scattered (Sedna and (26181)), Resonant 
((47171), (55638), (84922) and (208996)) and Classical ((19308) and Quaoar) bodies. Thus, 
our survey and previous spectroscopic works independently suggest that ices are found in 
among all orbital classes throughout the Kuiper Belt and are thus not uniformly destroyed 
by solar UV photons and cosmic rays. Nor does it appear that ice surfaces are a extreme 
minority exception to the radiation processing scenario, as roughly 10% of the non-Haumea 
bodies in our survey have ice of some kind on their surfaces. 



5.2. Sensitivity 

Our sensitivity to ices can be estimated by comparing the results of our photometry 
with spectroscopic observations from literature, as seen in Table |3l For all of the Haumea 
family members, which exhibit extremely large water ice absorption from 55% - 80% depth 
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in absorption, our photometric work and spectroscopic works agree. For non-Haumea bodies, 
some of which show very low levels of absorption (~ 15% absorption depth), the results are 
somewhat discrepant, which marks the minimum absorption levels to which we are sensitive. 
In addition, bodies with significant quantities of ultra-red matter may have ices that are not 
easily detectable with our J-band continuum estimation. 



5.2.1. Ultra-Red Matter and J-band Continuum Estimation 

A key quantity for estimating the amount of absorption in a spectrum is an accurate 
assessment of the continuum, i.e the reflectance found outside of any absorptions. For our 
survey, we have chosen the J-band because it is generally unaffected b y ices, as des cribed in 



§2.11 Unfortunately, the presence of the so-called "ultra- red matter" (jJewittll2002l ) can still 
impact this estimate. The ultra-red matter, which primarily colors the visible, does extend 
partway into the J-b and to about 1.3 ra icrons. The best example of this spectroscopically 



is found for Quaoar (jJewitt fc Luull2004l ). which is one of the reddest KBOs in the visible 
(and thus has a large amount of ultra-red matter), but has a relatively colorless near-IR 
continuum (attested by the fact that the near-ir shape can be fit with a simple water ice 
model with no broad coloring agent). A careful examination of Quaoar's spectrum indicates 
that the ultra-red matter likely extends into the J-band, although the exact point where its 
absorption ends is not clear, since the spectral shape of the ultra-red matter is unknown. 
However, if we assume that Quaoar's J-band spectrum differs from that of pure water ice 
only due to the contribution of ultra-red matter, then the maximum contamination that the 
ultra-red matter can impart to a J-band continuum measurement is 0.15 magnitudes, which 
is significantly smaller than the ~ 1 magnitude of color difference between neutral and icy 
bodies in our experiement. 

Thus, for Quaoar, and any other bodies with extremely red visible slopes, the J-band 
continuum measurement could be depressed by 0.15 magnitudes using our method. Thus, 
we use the Y — J colors in our survey to identify extremely red KBOs, which could harbor 
more ices on their surfaces than our photometry indicates. We have flagged these bodies 
in Table [3] by examining the F — J colors we have collected. We indicate red bodies as 
being those bodies with F — J > 0.1 to better than 3a confidence. Of the 29 bodies with 
Y — J known to better than 0.15 magnitudes, only 3 bodies have such red colors: Quaoar, 
2005 UQ513 and 2004 PF115. Thus, the point of rejecting objects with extreme Y — J values 
is that red Y — J color is a sign that the J-band flux measurement may be contaminated 
by ultra-red matter, and therefore ice detection may be depressed by not more than 0.15 
magnitudes. 
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5.2.2. Absorption Depth and Surface Ice Fraction 



Rejecting the reddest objects in the sample yields a consistent picture between spec- 
troscopy and photometry — our survey can identify objects with > 15% water ice absorption 
depth. Ixion and Orcus were the two non-Haumea family members identified as having wa- 
ter ice by our photometric survey with band depths also measured by spectr oscopy. They 
have ~ 1 0% and ~ 20% wate r ice absorption depth s res pectively as found b y iTrujillo et al. 



(120051), Ide Bergh et all (120051 ). iBarucci et al.l ( 120081 ) and iMerhn et all (120101 ). Seven bodies 
were suspected to have water ice in spectroscopic works, but didn't indicate water in our 
work, as shown in Table [31 The mean water ice band depth for these bodies is 12.5% ±2.5%. 
Thus, we estimate that for our mean sample color accuracy of 0.1 magnitudes, we are sen- 
sitive to water ice band depths > 15%. Since the bodies with large methane absorption 
show a similar ~ 1 magnitude separation from solar bodies in our system, we expect similar 
sensitivities for the methane band. 

We also estimated our sensitivity by producing a very simple model of a range of bodies 
from neutral, to pure water ice, to pure methane ice, as s hown in Figure HI We use a simple 
bi directional refl ectance model similar to that described in lTruiillo et al.l (120051. l2007a[). b ased 



on 



Hapkd (ll993f). Using absorpti on coefficients for water ice 



Grundy fc Schmittl ( 



19981 ) and 



methane ice iGrundy et al.l (120021 ) . we computed the full-disk albedo expected from a body 
covered in 1 mm diameter grains at 40 K. We computed a suite of models with these param- 
eters, as presented as the grey triangular matrix in Figure lH The vertical scale bar labeled 
"Ice Fraction" represents total ice percentage (methane ice and water ice combined), with 
values less than 100% indicating the pres ence of a fictitious neu tral absorber with albedo 
equal to the mean albedo of KBOs, 0.1 (IStansberry et al.l 120081 ). The diagonal scale bar 
labeled "Methane Fraction" indicates the percentage of methane ice, with values less than 
100% including a linear combination of water ice. Thus, the top corner of the triangular grid 
represents a neutral body, the lower rightmost corner represents a body of 100% methane 
ice, and the leftmost corner represents a body of 100% water ice. Although many of the 
bodies in the graph appear to fall outside of the grid, this is expected for our crude three 
component model, as the real variety of surfaces on KBOs is certainly much more complex. 
The model does have utility, however, in that it provides a rough constraint on the minimum 
surface fraction of ice to which we are sensitive, roughly ~ 10% surface fraction in our model, 
which is consistent with the ~ 15% water band depth for both Orcus and Ixion. 
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6. Discussion 

6.1. Size versus Surface Type 

We see two clear trends in our data set: (1) the surfaces of the smallest Haumea family 
members have a much stronger incidence of water ice than Haumea itself and (2) only the 
largest KBOs have methane ice. 



6.1.1. Methane on Large Bodies 



To date, there are only three bodies known in the Kuiper Belt with strong methane 
detections, Makemake and Eris, as shown by our survey in Figure |5] and Pluto. We use 
if as a proxy for size, which implicitly assumes a fixed albedo among the bodies, because 
the vast majority of bodies in the plot have unknown albedos. The pos s ibility that only 
large bodies may harbor methane ice was suggested by ISchaller fc Brownl (j2007bl ) and was 
used as the basis for a simple model explaining surface phenomena based on vapor pressure. 
Simply put, only very large bodies have surface gravity strong enough to overcome methane 
vapor pressure even at large heliocentric distances. We tested this theory measuring a large 
number of bodies with a variety of sizes with a single methodology. We find no strong 
evidence for methane on any bodies in our sample besides Makemake and Eris. Although 
both Quaoar and Sedna have been reported to have alkanes on their surfaces, neither has 



a large enough am o unt of the material to be detected in our suryev (iJewitt fc Luu 



2004 



Barucci et al.ll2005l : ISchaller fc Brownl l2007al : iTruiillo et al.ll2007bl : lEmerv et al.l 120071 ). For 

our sample of 51 objects, the 2 with the smallest absolute magnitudes (a proxy for largest 
size in the fixed albedo case) have methane ice. The probability that this could occur in 
size-neutral process which coated 2 of 51 bodies with methane ice is given by the binomial 
theorem as 0.12%, thus this possibility c an be rejected with about 3.5cr confidence. To date, 
our observations are consistent with the ISchaller fc Brownl (l2007bl ) theory. 



6.1.2. Water in the Haumea Family 

Examining water ice on the Haumea family members we find that water ice presence is 
inversely correlated with size, or the smallest bodies are the ones with the most water ice on 
their surface, as shown in Figure |6l In this figure, we use absolute magnitude, as a proxy 
for size since the albedos of Haumea's progeny are unknown. Using if as a proxy for size 
implicitly assumes a fixed albedo among the bodies, which may not be the case. However, 
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due to the large magnitude difference (~ 3 magnitudes, or a factor 16 in flux) between 
Haumea and its progeny, it is extremely unlikey that even the brightest of the fainter bodies 
are larger than Haumea. We do not believe this to be the result of any observational bias, as 
our standard star observations show no such correlation with apparent magnitude. We assess 
the statistical validity of this result by showing that the Haumea J — H2O = —0.85 ± 0.04 
with H = 0.2 is significantly different from the rest of the Haumea family, all of which have 
3.3 < H < 5.3 and a population mean of J — H2O = —1.5 ±0.1 (excluding 2009 YE7). Thus, 
Haumea has less water absorption than the rest of the Haumea family with a confidence 
level of ~ 6(7. 

It is not clear the underlying physical cause of such a difference, however, there are 
two obvious possibilities: (1) the smaller bodies have a larger fraction of water ice on their 
surfaces or (2) the smaller bodies have a larger grain size, leading to stronger absorption 
features. Both of these possibilities are physically plausible. If we assume that the Haumea 
family was created in a collisional event involv ing a much larger KB O, it is very likely 
that this body had some form of differentiation ( McKinnon et allboosl ). Thus, the largest 
remaining member, Haumea, likely contains a larger fraction of rock in its core than the 
smaller family members which were created from the ejecta material, which was presumably 
closer to the surface, less dense, and with less structural strength than the Haumea material. 
Thus it is easily plausible that scenario (1) is valid because the family members contain more 
water ice on a bulk level. Scenario (2) is also plausible since any ice transport is likely very 
strongly correlated with surface gravity. If significant surface evolution has taken place on 
the Haumea family, which has been s uggested as a possib le explanation for the crystalline 



water ice presence on their surfaces (INewman et al. 



20081 ). then it is very likely that this 



process is affected by surface gravity. Surface gravity profoundly affects ice transport and 
geologic processes such as cryovolcanism, thus nearly any evolutional p rocess whether is b e 
external (e.g. micro-meteorite b ombardment leading to ice transport (jPorter et al.ll2010l )) 
or internal (e.g. cryovolcanism (jjewitt fc Luull2004l )) in origin likely is affected by body 
size. It should be noted that although cryovolcanism may be expected for larger bodies, 
in-situ observations of the ~ 2700 km diameter Triton by Voy ager show only a small amount 
of surface is likely affected by cryovolcanism (jAbelsonlll989l ). Recent theoretical modeling 
shows that cryovolcanism in large (> 1200 km diameter) KBOs may be enhanced over that 
of Triton (jPesch et al.l |2009| ) , making the role of cryovolcanism in disk- integrated surface 
properties still an open question. 
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6.2. The Role of Photometry in Future Work 

Our custom photometric survey is not meant to replace spectroscopy in any way. The 
limited spectral information provided in a photometric survey covering only a few filters 
doesn't provide anywhere near the amount of information needed to disentangle the astro- 
physical processes taking place on minor bodies in the outer solar system. Indeed, in-situ 
observations of Saturnian satellites such as lapetus, Phoebe and Hyperion show how remark- 
ably diverse small body surfaces are in the outer solar system with extremely high signal- 
to-noise measurements. Thus, our survey is not meant to be an ending point of Kuiper Belt 
surface study but a starting point that will identify the bodies of greatest interest for later 
follow-up work. 

Photometry does offer several distinct advantages over spectroscopy: 

1. Numerous instrumental losses are associated with spectroscopy. Although each instru- 
mental setup is different, observations from Gemini North telescope are fairly typical 
for near-infrared instruments. Losses associated with an hour of NIRI spectroscopy 
would typically include the following: 20% loss from light near the periphery of the 
object that is occulted by the slit, 10% acquisition overhead required to put the target 
in the slit, a 50% loss due to grism throughput and image quality degradation through 
the grism. In addition, spectral information is spread across a few thousand pixels of 
the detector, whereas photometric information is tightly associated in twenty pixels, 
greatly reducing instrumental noise. 

2. Sky line removal can be difficult for spectral integrations which often require ~ 5 
minute exposures for faint object signal to be apparent. Since OH lines emission can 
typically change on ~ 5 minute timescales, integrations longer than this can make sky 
line removal challenging. Since photometric integrations are typically of order ~ 1 
minute, pairwise subtraction of images efficiently removes most of the changing sky 
background. 

3. In crowded fields, potential contamination of bodies can be easily identified, as depicted 
in Figure [3 Spectroscopic observations typically only spend ~ 10 minutes imaging 
the science target during acquisition, with spectroscopic observations taking the vast 
majority of the telescope time. Thus, it is impossible to construct a deep-field image 
of the science target environment when performing spectroscopy, which, especially for 
objects near the galactic plane, can contain many contaminant background stars. 

Because of these advantages, we find that we can save roughly a factor ~ 3 in telescope 
time while still identifying the objects of key interest, such as the Haumea family members. 
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This gives us a large potential increase in sample size, and allows us to observe many more 
objects with a consistent instrumental setup in consistent conditions than an equivalent 
spectroscopic work. In addition, a population-wide study such as ours has the ability to 
identify important population outliers for further spectroscopic study by the sheer number 
of objects available in our sample. 

In this work, our primary goal is to present the details of the custom photometric filter 
system. A future work will include a must larger sample of bodies, including faint bodies 
that have to date been out of reach of spectroscopic works. Any bodies with unusual surface 
features will be identified for further study in a subsequent work. Since many of these bodies 
will be fainter than the V ^ 22 brightness level, they will require 10 ~ 30 hours of 10 meter 
telescope time for follow-up spectroscopy to determine more about their surfaces. Requesting 
this amount of telescope time for a single object with no prior information about its surface 
is simply not feasible. Nor is it a wise use of telescope resources when most faint bodies 
appear to be bland. Thus, we believe that future results from our survey will become the 
cornerstone of spectroscopic study of faint bodies in the upcoming decade and as 30 meter 
class telescopes begin to come on line. 

7. Summary 

The primary purpose of this work is to detail our custom photometric system and 
demonstrate that it is useful in identifying dynamical correlations on the subset of KBOs we 
have observed to date. We find the following in this work. 

1. Our custom photometric work can identify KBO surfaces with a factor ~ 3 less tele- 
scope time than similar quality spectroscopic works. We can discriminate neutral 
bodies from water and methane ice bodies with > 15% absorption depths, which cor- 
responds to approximately > 10% surface fractions of ices for simple models of KBO 
surfaces. 

2. We identify 9 Haumea family members including a new object 2009 YE7, which is 
consistent in visible color, water ice depth and proper orbital elements with the rest 
of the Haumea family. We reject several objects with dynamics similar to Haumea 
because they do not have water ice on their surfaces. 

3. We find that the smallest Haumea family members systematically have deeper water 
ice absorption than Haumea to 6a confidence. This observation is consistent with a 
coUisional origin for the system. 
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4. We identify two non-Haumea bodies with water ice in 3:2 mean-motion resonance 
with Neptune (Orcus and Ixion). Combining this with previous spectroscopic studies 
demonstrates that moderate water ice is pervasive throughout all Kuiper Belt dynam- 
ical classes. 

5. Outside the Haumea family, we find no evidence for any correlation between body size 
and water ice fraction. 

6. We find that only the largest KBOs harbor methane ice, which is consistent with 
ar guments of volatile loss tim escales based on surface gravity, such as those presented 



bv lSchaller fc Brownl (]2007b[ l. 
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Table 1. Filter Wavelength Definitions 



Filter 


-^min 


-^max 


-^central 


AA 








[fiUl] 


[urn] 


Water 


1.9490 


2.1195 


2.0343 


0.1705 


Methane 


2.1990 


2.3505 


2.2748 


0.1515 



Note. — Wavelengths are defined to be the 
wavelength where the filter transmissions reached 
50% of the maximum transmission. Minimum and 

maximum wavelengths are given by Amin and Amaxj 
respectively while central wavelength and filter 
widths are given by Acentrai and AA, respectively. 



Table 2. Observational Circumstances 



\JO}GCZ 


wDject 


Provisional 


oite 


u i uate 


UT Time 




Airmass 


J 


±12 u 


L^rl4 


Calibration 




Number 


Name 


Designation 






Range 




Range 






pj 


Sources 








2000 CNios 


Gemini 


2009-02-01 


12:25 - 


13:37 


1, 


.014 - 


1, 


.065 


450 


1440 


1320 


HIP24336, HIP55398 








2001 KD77 


Gemini 


2009-06-15 


12;16 - 


13:06 


1, 


,560 


1 


.887 


150 


360 


240 


HIP90869 








2001 QC298 


Gemini 


2009-08-04 


12;29 - 


13:27 


1, 


,049 - 


1 


.078 


180 


780 


420 


HIP111063, HIP117367, HIP9829 






2002 GH32 


Gemini 


2009-01-31 


15;30 - 


16:08 


1, 


.223 - 


1 


.259 


270 


660 


720 


HIP23259, HIP49580, HIP49942 






2002 KW14 


Gemini 


2009-06-30 


09:34 - 


10:34 


1, 


.515 - 


1, 


.887 


360 


1080 


1080 


HIP78107, HIP80609 








2002 MS4 


Magellan 


2009-04-05 


09:09 - 


10:34 


1, 


.066 - 


1, 


.113 


525 


600 


1170 


HD110747, HD118928, HD138159 






























HD142801, HD170717, HD95868 






2002 XV93 


Gemini 


2009-01-31 


10:05 - 


10:43 


1, 


.455 - 


1, 


.715 


270 


660 


720 


HIP23259, HIP49580, HIP49942 






2003 FE128 


Gemini 


2009-04-15 


12:33 - 


13:13 


1, 


.241 - 


1 


.336 


270 


720 


600 


HIP73606 








2003 UZ117 


Gemini 


2009-09-20 


13:13 - 


14:11 


1, 


.032 - 


1 


.049 


360 


540 


960 


HIP 19767, HIP7373 








2003 UZ117 


Gemini 


2009-12-30 


06:11 - 


07:07 


1, 


.035 - 


1 


.071 


360 


960 


1020 


HIP117367, HIP18768, HIP19767 






2003 UZ413 


Gemini 


2009-01-09 


07:18 - 


08:01 


1, 


.105 - 


1, 


.205 


270 


660 


960 


HIP 19767 








2004 NT33 


Gemini 


2009-06-24 


12:55 - 


13:50 


1, 


.004 - 


1, 


.025 


360 


1080 


1080 


HIP106356 








2004 NT33 


Gemini 


2009-08-06 


10:39 - 


11:36 


1, 


.012 - 


1 


.074 


330 


960 


1080 


HIP106356 








2004 PT107 


Gemini 


2009-09-10 


06:00 - 


06:37 


1, 


.411 - 


1 


.603 


270 


660 


540 


HIP111063, HIP117367 






2005 CB79 


Gemini 


2009-01-31 


11:47 - 


12:26 


1, 


.019 - 


1 


.054 


270 


600 


660 


HIP23259, HIP49580, HIP49942 






2005 GEi87 


Gemini 


2009-01-31 


13:58 - 


14:36 


1, 


.427 - 


1, 


.655 


270 


720 


420 


HIP23259, HIP49580, HIP49942 






2005 QU182 


Gemini 


2009-09-20 


08:40 - 


09:57 


1, 


.218 - 


1, 


.499 


390 


1440 


1380 


HIP19767, HIP7373 








2007 JH43 


Magellan 


2009-04-04 


07:33 - 


09:19 


1, 


.022 - 


1, 


.078 


600 


1320 


1935 


HD123385, HD138159, 


HD154805 






























HD77533, HD85538 








2008 LP17 


Gemini 


2009-06-29 


07:39 - 


08:51 


1 


.183 - 


1 


.361 


450 


1440 


1380 


HIP113050, HIP75923 








2009 YE7 


Gemini 


2009-12-30 


07:20 - 


07:59 


1, 


.146 - 


1 


.154 


270 


480 


660 


HIP117367, HIP18768, 


HIP19767 






2009 YE7 


Gemini 


2010-02-21 


05:42 - 


06:35 


1, 


.231 - 


1 


.405 


330 


720 


420 


HIP18768 




(19308) 




1996 TOee 


Gemini 


2009-09-18 


09:56 - 


10:45 


1, 


.011 - 


1, 


.048 


270 


660 


660 


HIP117367, HIP11747, 


HIP9829 


(19308) 




1996 TOee 


Gemini 


2009-12-30 


05:23 - 


06:03 


1, 


.058 - 


1, 


.132 


270 


600 


540 


HIP117367, HIP18768, 


HIP19767 


(24835) 




1995 SMss 


Gemini 


2009-09-18 


13:13 - 


14:07 


1, 


.007 - 


1, 


.044 


270 


1260 


660 


HIP117367, HIP11747, 


HIP9829 


(26181) 




1996 GQ21 


Gemini 


2009-06-30 


07:49 - 


09:01 


1, 


.193 - 


1, 


.395 


450 


1440 


1260 


HIP78107, HIP80609 




(26375) 




1999 DEg 


Magellan 


2009-04-05 


02:09 - 


04:07 


1, 


.140 - 


1, 


.223 


450 


1980 


1665 


HD110747, HD118928, 


HD138159 






























HD142801, HD170717, 


HD95868 


(28978) 


Ixion 


2001 KX76 


Gemini 


2009-05-01 


14:43 - 


15:03 


1, 


.664 - 


1 


.790 


180 


300 


300 


HIP83875, HIP84181 




(38628) 


Huya 


2000 EB173 


Magellan 


2009-04-05 


05:06 - 


05:59 


1, 


.152 - 


1, 


.263 


300 


600 


990 


HD110747, HD118928, 


HD138159 






























HD142801, HD170717, 


HD95868 


(40314) 




1999 KR16 


Gemini 


2009-01-31 


14:47 - 


15:26 


1, 


.185 - 


1 


.263 


270 


720 


660 


HIP23259, HIP49580, HIP49942 


(47171) 




1999 TC36 


Magellan 


2008-10-18 


04:03 - 


05:03 


1, 


.129 - 


1, 


.148 


450 


810 


840 


HD9729, HIP11747 





Table 2 — Continued 







Provisional 


oite 


UT Date 


UT Time 


Airmass 


T 

J 


±12 w 


Orl4 


Oalibration 


Number 


Name 


Designation 








Range 


Range 








Sources 


(47932) 




2000 


GNi7i 


Gemini 


2009-05- 


■16 


10:44 


- 11:24 


1.384 - 


1 


.550 


270 


720 


660 


HIP73606 


(50000) 


QucLoax 


2002 


LMeo 


Magellan 


2009-04- 


■05 


07:23 


- 07:54 


1.111 - 


1 


.181 


300 


300 


300 


HD110747, HD118928, HD138159 
































HD142801, HD170717, HD95868 


(55565) 




2002 


AWl97 


Gemini 


2009-01- 


■31 


12:40 


- 13:18 


1.146 - 


1 


.261 


270 


660 


720 


HIP23259, HIP49580, HIP49942 


(55565) 




2002 


AWl97 


Magellan 


2009-04- 


■04 


02:06 


- 02:57 


1.229 - 


1, 


.312 


450 


660 


660 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(55636) 




2002 


TX300 


Gemini 


2008-11- 


■06 


06:49 


- 07:21 


1.015 - 


1 


.037 


180 


720 


300 


HIP394 


(55636) 




2002 


TX300 


Gemini 


2009-01- 


■10 


04:35 


- 04:53 


1.041 - 


1 


.065 


360 


300 


300 


HIP394 


(55637) 




2002 


UX25 


Magellan 


2008-10- 


■18 


05:24 


- 06:24 


1.286 - 


1 


.366 


450 


540 


1080 


HD9729, HIP11747 


(55638) 




2002 


VE95 


Gemini 


2009-02- 


■01 


07:36 


- 08:31 


1.127 - 


1 


.304 


360 


1080 


1080 


HIP24336, HIP55398 


(65489) 


Ceto 


2003 


FX128 


Gemini 


2009-05- 


■13 


08:00 


- 08:52 


1.061 - 


1 


.129 


240 


780 


660 


HIP72855, HIP82233 


(84522) 




2002 


TC302 


Gemini 


2009-08- 


■04 


13:55 


- 14:36 


1.028 - 


1 


.086 


270 


720 


720 


HIP111063, HIP117367, HIP9829 


(84922) 




2003 


VS2 


Gemini 


2008-12- 


■02 


10:49 


- 11:27 


1.081 - 


1 


.147 


270 


720 


720 


HIP21333 


(86177) 




1999 


RY215 


Gemini 


2009-09- 


■19 


11:08 


- 12:23 


1.225 - 


1 


.595 


420 


1260 


1260 


HIP117367 


(90377) 


Sedna 


2003 


VB12 


Gemini 


2008-12- 


■03 


09:54 


- 11:06 


1.060 - 


1 


.209 


450 


1440 


1440 


HIP18768 


(90482) 


Orcus 


2004 


DW 


Magellan 


2009-04- 


■04 


00:32 


- 01:40 


1.091 - 


1 


.159 


300 


1080 


1080 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(90568) 




2004 


GV9 


Magellan 


2009-04- 


■04 


06:10 


- 07:15 


1.001 - 


1 


.032 


900 


300 


480 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(119951) 




2002 


KX14 


Gemini 


2009-05- 


■01 


13:22 


- 14:09 


1.441 - 


1 


.629 


270 


720 


720 


HIP83875, HIP84181 


(120132) 




2003 


FY128 


Magellan 


2009-04- 


■04 


03:05 


- 04:11 


1.067 - 


1 


.163 


450 


660 


1260 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(120178) 




2003 


OP32 


Magellan 


2008-10- 


■18 


01:23 


- 03:42 


1.205 - 


1 


.727 


930 


1200 


2220 


HD9729, HIP 11747 


(120347) 




2004 


SB 60 


Gemini 


2008-11- 


■28 


04:53 


- 06:08 


1.006 - 


1, 


.085 


450 


1440 


1440 


HIP113688 


(120348) 




2004 


TY364 


Gemini 


2009-09- 


■18 


12:00 


- 12:40 


1.167 - 


1 


.196 


270 


720 


720 


HIP 117367, HIP 11747, HIP9829 


(136108) 


Haumea 


2003 


ELsi 


Magellan 


2009-04- 


■05 


04:24 


- 04:50 


1.574 - 


1 


.652 


150 


300 


450 


HD110747, HD118928, HD138159 
































HD142801, HD170717, HD95868 


(136108) 


Haumea 


2003 


ELai 


Magellan 


2009-04- 


■04 


05:22 


- 05:52 


1.505 - 


1 


.525 


300 


300 


300 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(136199) 


Eris 


2003 


UB313 


Magellan 


2008-10- 


■18 


06:30 


- 07:16 


1.215 - 


1 


.371 


300 


300 


1020 


HD9729, HIP11747 


(136472) 


Makemake 


2005 


FY9 


Magellan 


2009-04- 


■04 


04:34 


- 05:05 


1.881 - 


1, 


.891 


150 


300 


510 


HD123385, HD138159, HD154805 
































HD77533, HD85538 


(145451) 




2005 


RM43 


Magellan 


2008-10- 


■18 


07:36 


- 09:09 


1.213 - 


1 


.445 


750 


1260 


1380 


HD9729, HIP11747 


(145451) 




2005 


RM43 


Magellan 


2008-10- 


■19 


08:43 


- 09:32 


1.359 - 


1 


.584 


300 


810 


1080 


HIP 18768 



Table 2 — Continued 







Provisional 


oite 


u i uate 


UT Time 


Airmass 


T 
J 


-ri2 w 


C'-ri4 


Oali Drat ion 


Number 


Name 


Designation 






Range 




Range 






la] 


Sources 


(145452) 




2005 RN43 


Gemini 


2008-11-01 


05:22 - 


06 


36 


1.067 - 


1.116 


450 


1440 


1440 


HIPl 10035, HIP394 


(145453) 




2005 RR43 


Magellan 


2008-10-19 


04:15 - 


05 


58 


1.203 - 


1.498 


600 


1890 


1260 


HIP18768 


(145480) 




2005 TB190 


Gemini 


2009-08-07 


09:34 - 


10 


48 


1.119 - 


1.338 


420 


1440 


1380 


HIP117367, HIP80609 


(145480) 




2005 TB190 


Gemini 


2009-09-10 


07:12 - 


08 


26 


1.132 - 


1.371 


450 


1440 


1440 


HIP111063, HIP117367 


(174567) 




2003 MW12 


Magellan 


2009-04-05 


06:17 - 


07 


16 


1.235 - 


1.451 


450 


660 


990 


HD110747, HD118928, HD138159 
HD142801, HD170717, HD95868 


(175113) 




2004 PF116 


Gemini 


2009-08-08 


10:56 - 


12:12 


1.378 - 


1.407 


450 


1440 


1440 


HIP110512 


(202421) 




2005 UQ513 


Gemini 


2008-11-01 


07:04 - 


08:19 


1.009 - 


1.023 


450 


1440 


1440 


HIP110035, HIP394 


(208996) 




2003 AZ84 


Gemini 


2009-01-31 


11:01 - 


11:39 


1.124 - 


1.236 


270 


720 


720 


HIP23259, HIP49580, HIP49942 


(208996) 




2003 AZs4 


Magellan 


2009-04-05 


00:23 - 


01:53 


1.359 - 


1.630 


600 


1320 


1440 


HD110747, HD118928, HD138159 




























HD142801, HD170717, HD95868 


(225088) 




2007 ORio 


Gemini 


2009-06-29 


12:38 - 


13:50 


1.213 - 


1.321 


450 


1440 


1380 


HIPl 13050, HIP75923 


(225088) 




2007 ORio 


Gemini 


2009-09-16 


08:10 - 


09:25 


1.212 - 


1.239 


450 


1380 


1440 


HIP113050 


(229762) 




2007 UK126 


Magellan 


2008-10-19 


06:23 - 


08:23 


1.113 - 


1.184 


360 


2250 


1380 


HIP18768 



Table 3. Photometry 



Object 


Object 


Provisional 


J 


J - H2O 


J - CHi 


Photometric 


Spectroscopic 


Number 


Name 


Desif 


^nation 


[maerl 


fmaefl 


[maerl 


Tvne 


Absorption^ 






2000 


CNio5 


20.48 ± 0.06 


0.09 ± 0.09 


-0.28 ±0.20 


Neutral 








2002 


KWi4 


19.69 ±0.06 


-0.05 ±0.09 


-0.16±0.17 


Neutral 








2002 


MS4 


19.71 ±0.27 


0.18 ±0.33 


0.11 ±0.30 


Neutral 








2002 


XV93 


19.80 ±0.17 


0.13 ±0.19 


0.07 ±0.34 


Neutral 








2003 


FE128 


20.54 ±0.13 


0.19 ±0.16 


0.09 ±0.35 


Neutral 








2003 


UZny 


20.27 ±0.07 


-1.25 ±0.27 


-0.96 ±0.23 


Water^ 


~ 80% Water^ 






2003 


UZ413 


19.59 ±0.07 


0.08 ±0.11 


-0.06 ±0.12 


Neutral 








2004 


NT33 


19.17 ±0.04 


-0.11 ±0.08 


-0.00 ±0.09 


Neutral 








2004 


PT107 


20.23 ±0.11 


0.11 ±0.17 


0.07 ±0.23 


Neutral 








2005 


CB79 


19.61 ± 0.04 


-1.84 ±0.25 


-0.70 ±0.14 


Waterf 


~ 70% Water2 






2005 


QU182 


19.66 ±0.28 


0.10 ±0.29 


0.14 ±0.29 


Neutral 








2007 


JH43 


19.36 ± 0.04 


0.08 ±0.07 


-0.18 ±0.10 


Neutral 








2008 


LP17 


19.74 ±0.06 


-0.03 ±0.11 


0.01 ±0.10 


Neutral 








2009 


YE7 


20.84 ±0.11 


-1.03 ±0.27 


-0.63 ±0.40 


Water 




(19308) 




1996 


TOee 


20.73 ±0.23 


-1.20 ±0.40 


-0.06 ±0.29 


Water^ 


Wateri2.i4 


(24835) 




1995 


SM55 


19.43 ±0.05 


-1.87 ±0.26 


-0.88 ±0.18 


Water^ 


56% Wateri 


(26181) 




1996 


GQ21 


19.12 ±0.04 


0.23 ±0.05 


0.03 ±0.07 


Neutral 


9% Wateri 


(26375) 




1999 


DEg 


19.25 ±0.05 


-0.02 ±0.08 


-0.05 ±0.11 


Neutral 


Neutral^'^ 


(28978) 


Ixion 


2001 


KX76 


18.35 ±0.07 


-0.52 ±0.25 


0.01 ± 0.08 


Water 


~ 10% Water3'i3 


(38628) 


Huya 


2000 


EB173 


18.29 ±0.03 


0.06 ±0.06 


0.24 ±0.05 


Neutral 


Neutral^'^^ 


(40314) 




1999 


KR16 


19.48 ± 0.03 


0.01 ±0.04 


0.26 ±0.17 


Neutral 




(47171) 




1999 


TC36 


18.00 ±0.02 


0.07 ±0.03 


0.10 ±0.04 


Neutral 


8% Water! 


(47932) 




2000 


GN171 


19.07 ±0.04 


0.02 ±0.09 


0.09 ±0.08 


Neutral 


Neutral^ 


(50000) 


Quaoar 


2002 


LMeo 


17.25 ±0.02 


0.11 ±0.03 


0.27 ±0.04 


NeutraF 


22% Water4^6,io 


(55565) 




2002 


AWig7 


18.64 ±0.04 


0.09 ±0.05 


0.18 ±0.06 


Neutral 


Neutral^'^ 


(55636) 




2002 


TX300 


18.89 ±0.05 


-1.38 ±0.22 


-0.53 ±0.19 


Water* 


64% Water! 


(55637) 




2002 


UX25 


18.48 ±0.05 


-0.02 ±0.06 


0.12 ±0.05 


Neutral 


Neutral^ 



Table 3 — Continued 



Object 


Object 


Provisional 


J 


J-H2O 


J-CH4 


Photometric 


Spectroscopic 


Number 


Name 


Designation 


[mag] 


[mag] 


[mag] 


Type 


Absorption'^ 


(55638) 




2002 VE95 


18.22 ±0.02 


-0.08 ±0.03 


-0.06 ±0.03 


Neutral 


9% Wateri 


(65489) 


Ceto 


2003 FX128 


19.61 ±0.06 


0.09 ±0.08 


0.23 ±0.13 


Neutral 


14% Wateri 


(84522) 




2002 TC302 


18.93 ±0.02 


-0.02 ±0.06 


-0.13 ±0.05 


Neutral 


Neutral^ 


(84922) 




2003 VS2 


18.19 ±0.03 


0.20 ±0.05 


0.11 ±0.08 


Neutral 


7% Wateri 


(90377) 


Sedna 


2003 VB12 


19.06 ±0.04 


0.02 ±0.05 


-0.17 ±0.09 


Neutral 


~ 10% Water'' 


(90482) 


Orcus 


2004 DW 


18.12 ±0.03 


-0.51 ±0.04 


-0.13 ±0.05 


Water 


~ 20% Water^ 


(90568) 




2004 GVg 


18.55 ±0.03 


0.14 ±0.06 


0.29 ±0.07 


Neutral 




(119951) 




2002 KX14 


19.04 ±0.07 


-0.04 ±0.09 


0.13 ±0.10 


Neutral 


Neutral"'^ 


(120132) 




2003 FY128 


19.09 ±0.03 


0.04 ±0.09 


0.12 ±0.08 


Neutral 


Neutrali'3 


(120178) 




2003 OP32 


19.37 ±0.02 


-1.40 ±0.28 


-0.55 ±0.14 


Waterf 


74% Wateri 


(120347) 




2004 SBeo 


19.53 ±0.05 


0.22 ±0.06 


0.20 ±0.07 


Neutral 


NeutraP 


(120348) 




2004 TY364 


18.75 ±0.03 


-0.11 ±0.06 


-0.21 ±0.07 


Neutral 


Neutral -'^ 


(136108) 


Haumca 


2003 ELei 


16.73 ±0.03 


-0.85 ±0.04 


-0.25 ±0.05 


Watca-f 


- 60% Watcr"^ 


(136199) 


Eris 


2003 UB313 


17.99 ±0.03 


-0.01 ±0.05 


-1.57 ±0.12 


Methane 


~ 95% Methane^ 


(136472) 


Makemake 


2005 FYg 


16.56 ±0.02 


-0.18 ±0.04 


-1.59 ±0.07 


Methane 


~ 85% Methane^ 


(145451) 




2005 KM43 


19.08 ± 0.02 


—0.22 ± 0.06 


A AC 1 A AA 

0.05 ± 0.09 


iN eutral 




(145452) 




2005 RN43 


18.43 ±0.02 


0.11 ±0.03 


0.08 ±0.03 


Neutral 




(145453) 




2005 RR43 


19.20 ±0.02 


-1.38 ±0.17 


-0.62 ±0.15 


Water* 


65% Wateri 


(145480) 




2005 TB190 


19.44 ±0.04 


0.00 ±0.07 


-0.08 ± 0.09 


Neutral 




(174567) 




2003 MW12 


18.89 ±0.02 


0.08 ±0.06 


-0.08 ±0.07 


Neutral 




(175113) 




2004 PF115 


18.71 ±0.02 


0.02 ± 0.03 


-0.04 ±0.04 


NeutraF 




(202421) 




2005 UQ513 


19.04 ±0.04 


0.10 ±0.05 


0.06 ±0.06 


Neutral"" 




(208996) 




2003 AZ84 


19.37 ±0.03 


-0.11 ±0.08 


0.08 ±0.11 


Neutral 


18%,~ 25% Wateri'3 


(225088) 




2007 ORio 


19.06 ±0.03 


-0.15 ±0.05 


-0.08 ±0.06 


Neutral 





^Reported absorption refers to the percentage depth of the absorption feature, since surface fraction is a highly model- 
dependent quantity. 



^Likely Haumea family member based on composition (this work) and dynamics ( Ragozzine &: Brownll2007 ). 



"^We find these three objects to be extremely red in the visible, indicating that a continuum measurement is difficult to 
obtain. These objects may have more absorption than indicated from the J — H2O or J — CH4 metrics. We flag objects as 
"red" if they have F — J > 0.1 to at least Scr significance as determined from our photometry. See Section [5.2.1l for details. 



Note. — This table lists all objects observed and their surface types based upon our custom photometric system (this 
work). Spectroscopic surface types are based on other works as referenced in the table. Haumea family me mbers are iden 
tified as be ing both water-ice rich a nd ha ving orbital p arame ters s imilar to Haumea. References — (1) 



( 2nn8h. ('21 ISchaller fc Brownl (I2nn8h. (3) iMerlin et al.l (l20inh. m iTruiiUo et aP (l2007al). (5) 



Dalle Ore et all hoOflh. (7 ) 



Truiillo et al. 



2007 



(81 



Brown et al 



(I20n5h. (9) iTruiillo et all toO^. (10) | ,Tewitt fc Lmi 



Barkume et al 



Brown et al. (l2007al). (6) 



tooi ). riD l.Tewitt fc Lm3 |200ll ). (;i2l lBrown et all (|2007bll . (U) \de Bergh etall (|2005h . fM^ lNoh et al.l (j2000^ 
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Table 4. Ice Detections From This Work 



Object 


Object 


Provisional 


a* 


e* 


i* 


Dynamical 


Number 


Name 


Designation 


[AU] 




n 


Class 




Haumea family 






2003 UZii7 


44.26 


0.13 


27.88 


Haumea 






2005 CB79 


43.27 


0.13 


27.17 


Haumea 


(19308) 




1996 TOge 


43.32 


0.12 


28.02 


Haumea 


(24835) 




1995 SM55 


41.84 


0.10 


26.85 


Haumea 


(55636) 




2002 TX300 


43.29 


0.13 


26.98 


Haumea 


(120178) 




2003 OP32 


43.24 


0.10 


27.05 


Haumea 


(136108) 


Haumea 


2003 ELei 


43.10 


0.19 


26.85 


Haumea 


(145453) 




2005 RR43 


43.27 


0.13 


27.07 


Haumea 






2009 YE7 


44.43 


0.13 


27.99 


Haumea 


Non-Haumea Water Ice 


(90482) 


Orcus 


2004 DW 


39.17 


0.23 


20.56 


3:2 Resonance 


(28978) 


Ixion 


2001 KX76 


39.62 


0.25 


19.65 


3:2 Resonance 


Methane objects 


(136199) 


Eris 


2003 UB313 


68.01 


0.43 


43.83 


Scattered 


(136472) 


Makemake 


2005 FYg 


45.36 


0.16 


29.01 


Classical 



*For the Haumea family members except for 2009 YE7, a, e and i are th e proper orbital ele- 
ments semimajor axis, eccentricity and inclination, respectively, as presented bv lRagozzine fc Brown 
(120071 For 2009 YE7, proper orbital elements are presented for the first time in this work, and are 
an average of 20 million years of orbital motion computed by the MERCURY symplectic integrator 
(IChamberg|ll999h . The orbital elements for the non-Haumea family members are osculating elements 
provided by the Minor Planet Center. 



Note. — Bodies from our survey with significant amount of ices on their surfaces. The top section 
of the table are all Haumea family members based on their (1) wate r ice absorption and (2) th e 
proximity of their proper orbital elements to Haumea as identified by iRagozzine fc Brownl (|2007l l , 
with the exception of 2009 YE7 which is identified in this work based on the proper orbital elements 
we computed, the water ice absorption depth and the optical colors of the body. Orcus and Ixion 
are two non-family members with significant water ice absorption. 
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Table 5. Near- infrared Observations of KBOs Near the 
Haumea Family 



Object 


Object 


Provisional 






Water 


Number 


Name 


Designation 


[m s^-^] 




Ice 


(136108) 


Haumea 


2003 ELgi 





2.83 


Yes 


(19308) 




1996 TOee 


15.0 


2.83 


Yes 


(145453) 




2005 RR43 


58.0 


2.84 


Yes 






2003 UZ117 


60.8 


2.84 


Yes 






2005 CB79 


66.5 


2.84 


Yes 


(55636) 




2002 TX300 


68.4 


2.84 


Yes 


(120178) 




2003 OP32 


91.4 


2.85 


Yes 






2009 YE7 


- 100 


2.84 


Yes 


(24835) 




1995 SM55 


123.3 


2.84 


Yes 


(202421) 




2005 UQ513 


39.0 


2.84 


No 


(136472) 


Makemake 


2005 FYg 


118.0 


2.84 


No 






2004 PT107 


161.9 


2.83 


No 


(120347) 




2004 SBeo 


218.5 


2.86 


No 


(40314) 




1999 KR16 


242.9 


2.84 


No 



*As computed in iRagozzine fc Brownl (|2007[ ). excepting 2009 YE7. The 
explicit computation of Afmin for 2009 YE7 is beyond the scope of 
this work. However we include a crude estimate based on the fact 
that the proper orbital elements of 2009 YE7 are close to those of 
2003 UZ117 and wel l withi n the Ai/ = 150 m s~^ region of Figure 1 from 
Ragozzine fc Brownl |2007h . We find in this work that 2009 YE7 is a 
Haumea family member. 



Note. 



We list the bodies we observed that were also identified in 



Ragozzine fc Brownl (|2007l) as potential family members based on their dy- 
namical association with Haumea, ordered by Ai^niin- The bodies with 
water ice detected appear to have lower Afmin compared to those that 
have no water ice, although this result is not significant at the 3ct level. 
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Table 6. Visible Photometry for 2009 YE7 From This 

Work 



Bandpass Magnitude 







21.55 ±0.03 


9'- 


r' 


0.45 ±0.04 




%' 


0.30 ±0.05 


R 




21.35 ±0.03* 


B - 


R 


1.03 ±0.04* 


V - 


R 


0.38 ±0.04* 


R- 


I 


0.51 ±0.05* 


rriRi 


1,1,0°) 


4.13 ±0.05* 



* Sloan colors were converted to the Johnson-Morgan- 
Cousins BVRI colo r system using tr ansfer equations from 



Smith et al.l (|2002[ ). as described in ISheppardI ( l2010l ). 



Note. — At the time of the observations, 2009 YE7 was 
at heliocentric distance 50.665 AU, geocentric distance 
51.153 AU, phase angle 0.979 degrees. rriRil, 1, 0°) refers 
to a hypothetical red magnitude at lAU heliocentric and 
geocentric distance and phase angle of 0°. 
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Wavelength [microns] 



Fig. 1. — Our custom filter bands (grey regions) witli synthetic water and methane ice spectra 
overlaid (bottom two traces, offset vertically for clarity). Our primary goal is the discrimination 
between water ice, methane ice, and featureless KBO surfaces. The top two traces show sky 
absorption and sky OH emission, the two main sources of systematic error in the near-infrared. 
The Y filter and J filter are used to sample the continuum by being as wide as possible while still 
minimizing ice and telluric effects. The "Water" filter limits are set by telluric absorption on the 
blue end and the overlapping methane absorption on the red end. The "Methane" filter limits are 
set by overlapping water absorption on the blue end and thermal emission on the red end. The 
methane filter also has some sensitivity to other alkanes which have C-H absorption in the methane 
region. The top portion of the plot shows traditional HK bandpasses, a poor diagnostic for ices. 
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Fig. 2. — Efficiency curves for our two custom filters as a function of wavelength for the 
water ice filter (left) and the methane ice filter (right). 
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Fig. 3. — Photometric color calibration for Gemini/NIRI. 
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Fig. 4. — Combined photometry for all objects. Data symbol position is computed from our 
photometric survey. Data symbol shape is determined from published spectroscopy (see Table [3] for 
complete list of references). Objects clearly fall into three groups: Methane (Eris and Makemake), 
Water (Haumea family members now including 2009 YE7, Ixion, and Orcus) and Neutral / low 
absorption (all others). In general there is a strong correlation between our photometric work 
and published spectroscopic works for those bodies observed by both. The triangular grey grid 
represents colors expected for synthetic objects with varying percentages of methane, water and 
total ice fraction. The vertical axis of the grid labelled "Ice Fraction" represents total ice fraction 
(water and methane combined) with the remainder being represented by a neutral absorber. The 
diagonal axis labelled "Methane Fraction" represents the total fraction of ice that is methane ice, 
with water ice being the remainder of the ice fraction. The model details are discussed in full in 
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Fig. 5. — Methane ICG clS £l function of size using the proxies J — CH^ for methane ice 
fraction and H (absolute magnitude) as a proxy for size. Only the largest bodies in our 
sample, Makemake and Eris show evidence for methane ice. The dotted lines represent the 
±3(7 hmits to the non-methane bodies based upon their J — CHi values. Hollow circles are 
Haumea and its family members. 
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Fig. 6. — Water ice as a function of size using the proxies J — H2O for water ice fraction 
and H (absolute magnitude) as a proxy for size. The small Haumea family members show 
much larger water ice absorption than Haumea itself. This could be explained by (1) a larger 
surface fraction covered in water ice or (2) larger grain sizes. The dotted lines represent the 
±3(7 limits to the non-family members based upon their J — H2O values. Orcus and Ixion 
fall outside of these 3a limits and we thus consider them to have significant amounts of ice 
on their surfaces. Hollow circles are Haumea and its family members. 
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Fig. 7. — An example of a rich field from Magellan appears above. The KBO 2002 MS4 is 
shown in a deep composite image of all frames taken over a ~ 2 hour period, digitally tracked 
at the non-sidereal rate of the object. A deep image was constructed of every science field in 
our sample to ensure no stellar contamination. Such a verification procedure is not possible 
for traditional spectroscopic measurements and background stars can easily unknowingly 
contaminate observations. 



